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The field of Generalized Parton Distribution (GPD) benefited from a wealth of exclusive 
reactions measurements since 2000. Extraction of GPDs from observables has begun with this 
first generation of experiments. In the short and mid-term future the programs of COMPASS-II 
and JLab at 12 GeV will enlarge the existing kinematic domain and open an era of high- 
precision measurements. Tools for phenomenological analysis are being developed at the same 



Qh 

time and will be mature to handle these forthcoming measurements. In the long-term future 
the EIC project will allow an ambitious spin program in the small to intermediate xg domain. 
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Oh 1 Introduction 

Qh 

The Deeply Virtual Compton Scattering (DVCS) process was early recognized as the cleanest 
way to access Generalised Parton Distributions (GPD) and has been so far a very active 
field of research. However a flavour separation of GPDs will require Deeply Virtual Meson 
Production (DVMP) measurements (see reviews [1^] and references therein). 

One of the goals of this program is imaging the nucleon in 3d. Present observables may 
^ already be fit, giving indications on the actual sensitivity of observables to GPDs and hints 

^ to elaborate robust extracting methods. 

The first section is a brief reminder. The second part outlines some results of the extraction 
of GPDs as of 2011. The last section mentions promising future prospects. 



2 DVCS at leading twist and leading order 

Four GPDs H, E, H and E describe DVCS at leading twist and leading order, but the cross 
section depend on the Compton Form Factors (CFF) H, £, % and £ (see ref. [5]). The GPDs 
depend on the generalised Bjorken variable [5] £ (or equivalently on the standard Bjorken 
variable Xg), the virtuality of the initial photon Q 2 and the square momentum transfer 
between initial and final protons t. The CFFs are convolutions along the variable x of GPDs 
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with known kernels. This integration kernel yields a real and an imaginary part to a CFF. 
They are related by fixed- f dispersion relations [6,7] involving a subtraction constant related 
to the D-term [8]. However the D-term is poorly known and most of DVCS measurements 
were made in the region £ < 0.5. Using dispersion relations thus rely on models and may 
introduce biases in the extraction of GPDs from DVCS data. For that reason the real and 
imaginary parts of CFFs are taken as independent in some fitting procedures. 

3 Status of GPD analysis 

3.1 Counting degrees of freedom 

The problem consists in extracting four functions H, E, H and E of three variables 2 x, £ 
and t for each quark flavour (u, d and s). The knowledge of these functions should provide 
us Parton Distribution Functions (PDF) and Form Factors (FF). A naive counting quickly 
shows that building a flexible yet robust GPD parametrization is very involved because 
the number of free parameters has to remain low to keep the fit tractable and the numerics 
under control. 

3.2 Local fits 

This approach is detailed in ref . [9] and assumes the independence of the real and imaginary 
parts of CFFs. The mains assumptions are the validity of the twist-2 analysis of existing 
DVCS measurements and a negligible contribution of Im£. Each kinematic bin (xg, t, Q 2 ) 
is taken independently of the others, and the seven values Re %, Im H, Re £ , Im £, Re H, 
Im H and Re £ are extracted simultaneously. Nothing prevents large fluctuations between 
two neighbouring kinematic bins. In the following we will refer to these fits as local fits. 

The model-dependence is almost as low as possible but the problem is often ill-posed 
because the harmonic structure of a given observable {e.g. a BSA) at fixed kinematic usually 
depends on less than 7 parameters. Thus it is mandatory to combine different observables 
on the same kinematic bins (for example BSAs and Target Spin Asymmetries (TSA) or Beam 
Charge asymmetries (BCA)) and to assume boundaries for the values the CFFs may take. 

3.3 Global fits 

In the spirit of the work done on PDFs and FFs, global fits require a physically motivated 
parametrization of GPDs and deal with all observables on all kinematic bins at once. 

The free coefficients entering the expressions for GPDs are determined either from PDFs 
and FFs or from DVCS data. Two such studies have been reported so far for DVCS. In the 

2 We omit the Q 2 -dependence which is governed by the QCD evolution equations. 
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first case [10] all DVCS measurements on an unpolarized target are included in the fit in a 
(mostly) H-dominant framework. Fixed-f dispersion relations are used as a key ingredient. 
In the second case [11] CLAS and HERMES measurements are used in the fit, and JLab Hall 
A data are computed as a prediction of the fitting Ansatz ; all four GPDs H, E, H and E are 
described. 

The extrapolation outside the physical region (for example the interesting limits t — > or 
£ — > 0) is possible but an assessment of the parametrization uncertainties is a complex task 
in itself. 

3.4 Hydrid : Local / global fits 

This hydrid fitting procedure [12] is a combination of the previous two methods and has 
been applied with the main assumption of H-dominance and twist-2 accuracy. It involves 
a parametrization with the minimal properties : smoothness, leading order Q 2 evolution 
and polynomial parametrization. Since the functional is otherwise arbitrary it is a posteriori 
validated by the quality of the fit. 

The model dependence is tested by a systematic comparison to local fits and an estimate of 
the systematic error induced by the H-dominance hypothesis. The good agreement of the 
two methods (local fits vs global fits) is a necessary consistency check of this approach. The 
treatment of systematic errors is however a weak point since there is no rigorous procedure 
to evaluate them and only estimates are given using joint models and fits. 

3.5 Neural networks 

Neural network fits had been successfully performed for PDFs but their use for GPD extrac- 
tion is quite recent. First results are described in [13] within the H-dominance assumption. 
It is worth noting that it is a new development in the field of GPD extraction. It is hoped it 
will allow better estimates of errors when extrapolating outside the fitting domain. 

3.6 Key results 

To our knowledge, the first extraction of a combination of GPDs was made in [14] and the 
progress have been continuous so far. The different methods all agree on the following key 
results : 

• Dominance of twist 2 and validity of a GPD analysis of DVCS data. This is no trivial 
result considering the low Q 2 range accessible in Hermes and JLab kinematics. 

• ImH is the best determined quantity with an overall fair agreement between all 
estimates. However large uncertainties subsist on Re "H. 
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• There are already some indications about the invalidity of the H-dominance hypothe- 
sis even with the experimental statistical accuracy achieved so far. 

3.7 A trivial universality check 

Since GPDs have been modelled or extracted from DVCS and DVMP measurements, it is 
natural to check if these objects are identical. Although far from a satisfactory treatment 
of the universality question, as a first step we can take GPDs obtained in DVMP studies, 
predict DVCS observables and compare them to existing measurements. The predictions of 
the Kroll - Goloskokov model (see [15] and references therein) are depicted in Fig. 1 (left). In 
spite of the fair agreement between this model and the extracted values of H, the resulting 
DVCS prediction tends to overshoot the data. This is reminiscent of the current version of 
theVGG GPD model [1]. 
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Figure 1 : Measured CLAS BSAs for -0.4 < t < -0.2 GeV 2 . The red lines correspond to 
the prediction of the Kroll - Goloskokov GPD model (left). Projected CLAS12 BSAs for 
-0.35 < t < -0.25 GeV 2 . Statistical errors are smaller than data point bullets. The red 
lines correspond to a fit of the data, (right) 



4 Prospects 

4.1 Facilities and projects 

In 2012 COMPASS-II will start running its GPD program, giving access to several observ- 
ables with beam spin and charge differences, in between the kinematic domain of collider 
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and fixed-target experiments. 

By 2015 JLab's 12 GeV upgrade will reach a time of high-precision measurements with 
an expected few percent accuracy. See Fig. 1 (right) for projected data and a tentative fit. 
Despite an apparent agreement, the high accuracy of the data does not allow a fit assuming 
H-dominance any more. 

The science case for an electron-ion collider (EIC) has been developed since 2010 [16]. In the 
2020 decade we expect to produce a wealth of GPD related spin observables (see projected 
data in Fig. 2) with the aim of imaging the sea quark and gluon content of the nucleon. 
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Figure 2: Photon electroproduction beam spin asymmetries for the 20 x 250 EIC configu- 
ration, in the typical kinematic bin: 1.58 ■ 10~ 3 < x B < 2.51 ■ 10~ 3 , 3.16 < Q 2 < 5.61 GeV 2 
for four different f-bins as shown on each plot. Up to about 3 months of beam time with 
50% efficiency is necessary to achieve reasonable statistics in some of the bins. 



4.2 Phenomenological tools for GPD analysis 

A software platform is being developed to handle the phenomenological analysis of past and 
forthcoming data. It contains a database of experimental results, a database of theoretical 
predictions, a fitting engine and a visualizing software to compare measurements and 
model expectations. Modularity is a key-issue of the design : the same libraries are used to 
extract GPDs from measurements and to design new experiments connected to experimental 
descriptions and event generators. 

5 Conclusions 

During the 2000s GPD studies have made tremendous progress, with important measure- 
ments and encouraging first results on the extraction of GPDs. However a robust and 
efficient fitting strategy for DVCS and DVMP is still needed. 
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The 2010s will be very exciting times for GPDs : new facilities, high-precision measurements 
and bigger kinematic domains will constrain the phenomenology of GPDs. To this end a 
platform dedicated to global GPD analysis is currently under development. 
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